Optimized Dynamical Decoupling for Time Dependent Hamiltonians 
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The validity of optimized dynamical decoupling (DD) is extended to analytically time dependent 
Hamiltonians. As long as an expansion in time is possible the time dependence of the initial 
Hamiltonian does not affect the efficiency of optimized dynamical decoupling (UDD, Uhrig DD). 
This extension provides the analytic basis for (i) applying UDD to effective Hamiltonians in time 
dependent reference frames, for instance in the interaction picture of fast modes and for (ii) its 
application in hierarchical DD schemes with tt pulses about two perpendicular axes in spin space, 
to suppress general decoherence, i.e., longitudinal relaxation and dephasing. 

PACS numbers: 03.67.Pp, 82.56.Jn, 76.60.Lz, 03.65.Yz 



I. INTRODUCTION 

Progress in quantum information processing (QIP) re- 
quires a complete and coherent control of the dynamics of 
a quantum bit (spin S = 1/2) coupled to an environment 
(bath). In particular, one must be able to realize the 
no-operation reliably and coherently for long-time stor- 
age of quantum memory. Hence decoherence must be 
suppressed. The most general decoherence consists of 
both transversal dephasing and longitudinal relaxation, 
i.e., the decoherence rates l/'r2 1/Ti, respectively, 
in nuclear magnetic resonance (NMR) language. 

So far, only models without explicit time dependence 
have been considered to our knowledge. The terms in the 
Hamiltonian H without coherent control (we will call this 
Hamiltonian henceforth the initial one) do not have any 
explicit dependence on the time. For such a model tech- 
niques of various degrees of sophistication exist to sup- 
press the dephasing and/or the relaxationii^. We concen- 
trate here on the dynamical decoupling (DD)^"— which 
generalizes the original ideas on spin echo techniques to 
open systems and their application to QIP^^"^"'. Intu- 
itively, the interaction between spin (qubit) and bath is 
averaged to zero by means of repetitive tt pulses. Each 
pulse rotates the spin by an angle tt about a spin axis a, 
thus inverting its components perpendicular to a. 

A particularly efficient way to suppress pure dephas- 
ing is the optimized DD (Uhrig DD)2i"— li^"— where the 
instants tj (j G {1,2... N\) at which N instantaneous tt 



pulses^^ are applied are given by tj 
total time of the sequence and 



T^,- where T is the 



(5j = sin^(j7r/(2iV + 2)). 



(1) 



By efficient suppression it is meant that each pulse 
helps to suppress dephasing in one additional order in 
an expansion in T, i.e., N pulses reduce dephasing to 
0{T^^^\ The existence of an expansion in powers of 
T, at least as an asymptotic expansion, is a necessary 
assumption. 

So far, the derivation of the properties of UDD as de- 
fined in Eq. ([T|) was given for time independent initial 



Hamiltonians2r2£. The present study extends this deriva- 
tion to initial Hamiltonians including an analytic time 
dependence. This extension is a breakthrough because 
it establishes the applicability of optimized DD for ef- 
fective Hamiltonians in special reference frames, e.g., ro- 
tating frames, which induce an explicit time dependence. 
Such situations arise also where fast modes are treated in 
the interaction picture, are averaged over, or integrated 
out so that time dependent actions result and these ac- 
tions are sufficiently smooth in time. The condition on 
smoothness need not always be fulfilled. 

Another important application of UDD for time de- 
pendent Hamiltonians is the suppression of general deco- 
herence by the application of tt pulses around two per- 
pendicular spin axes on two hierarchical levels. If UDD 
worked only for time independent initial Hamiltonians 
the only known solution for the secondary level would 
be concatenation of primary UDD sequences^^. But re- 
cent numerical data by West et al. showed that also the 
suppression on the secondary level can be efficiently re- 
alized by UDD—. They called the scheme quadratic DD 
(QDD). Hence the derivation below provides the analytic 
foundation for the applicability of QDD. 

For the sake of simplicity, we first give the extended 
derivation for pure dephasing, addressing longitudinal re- 
laxation in a second step. Then the applications are dis- 
cussed again and we provide an explicit derivation of the 
time dependence of the effective Hamiltonian after the 
primary application of UDD. 



II. DEPHASING 

We consider the explicitly time dependent Hamiltonian 

H{t)^Hy,{t) + a,AAt) (2) 

is studied where the time dependences of the bath Hamil- 
tonian H\^{t) and of the coupling operator Az{t) are re- 
quired to be analytic, i.e., they can be expanded in t. 
If the system described by ^ is subject to N instan- 
taneous TT pulses at the instants {T5j},j e {1, 2, . . . , N} 
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about a spin axis perpendicular to the z-axis the effective 
Hamiltonian H{t) in the basis of unflipped spins reads 



H{t)^Hy,{t)+aMt)F{t) 



(3) 



where the switching function F{t) — ±1 appears which 
changes sign at the instants {T5j}. We are interested in 
the time evolution operator U (t) induced by H{t) 

U{T) := rexp(-i ( H{t)dt) (4a) 
= Uo{T)Ui{T) (4b) 

where T is the time-ordering operator. The second line 
is based on the interaction picture with respect to {t) 



Texp ( - i / Hh{t)dt) 
Jo 



Uo{t) 

Jo 

Ui{t) := rexp(-i / Ai{t)F{t)dt) 
Jo 



(5a) 

(5b) 
(5c) 



The key observation is that Ai{t) is analytic as well 
because H\j{t) is analytic according to our require- 
ment, and thus U(){t), and Az{t) again according to our 
requiremenliSl. Thus we have 



p=0 



(6) 



To be precise, to exclude any terms up to a given order 
iV we only need that Ai{t) can be represented by the 
sum 'yo. ^ wp to p — N plus a residual function of higher 
order, i.e., the convergence of the Taylor series is not 
needed. In general, the operators Ap are complicated 
integral expressions of the operators in ([2]). 

Next, the time evolution Ui{T) is expanded according 
to standard time dependent perturbation theory 

oo 

Un = / F{tn) ... F{h) / F{ti) 
Jo Jo Jo 

X Ai{tn)Ai{tn-l)...Ai{ti)dhdt2...dt^. (7b) 

Our aim is to show that the powers with n odd are of 
order T^+^ because only the odd powers in ct^ affect the 
qubit spin. Hence we can follow the reasoning of Yang 
and LiuF^ from here on. In order to keep the present 
communication self-contained we include the main steps. 
First, we expand in powers of T by inserting ^ into (I7bp 

Un ^ ^ T -^Pn ■ ■ ■ ^P2-^PlFpi,p2,...Pn (^) 

where pj G N and P„ :— J2'i=iPj ^^d 









^Pl...p 


„ := / din.. 


. f ' di2 f 




Jo 


Jo Jo 



We used the dimensionless relative times t := t/T. Since 
the N switching instants are given by TSj the function 
F{Tt) does not depend on T for given {5j}. Hence the 



coefficients K 



do not depend on T . 



Our goal is to show that Fp-^,,,p^ vanishes for n odd 
and N > ri + P„. Based on the UDD choice for the 
{Sj} in ([T|) the substitution i — sin^(6'/2) suggests itself 
because it renders f{9) :— F{T sir? (0 /2)) particularly 
simple if the {(5j} are chosen according to ((I]). Then 
/(6») = {~iy holds for 6* e {jir / {N + + 1)-k / {N +1) 
with j e {0, . . . , A^}. If we release this contraint on j 
allowing j € Z the function f{d) becomes an odd function 
with antiperiod tt/{N + I). Thus its Fourier series 

oo 

f{d) = E sin((2fc + l)(iV + 1)9) (10) 

k=0 

contains only harmonics sm{r9) with r an odd multiple 
of -I- 1. The precise coefficients C2k+i do not matter 
which can be exploited for other purposes, e.g., to deal 
with pulses of finite durationi^. 

Under the substitution i = sm^{9/2) the terms i^dt in 
© become sin^P(6'/2) sin(6')d6' which can be reexpressed 
as suitably weighted sum over terms sui{q9)d9 with q G 
^7 1 9 1 < P + 1- Thus we achieve our goal if we can show 
that the coefficients 

:= / d9n... [ ' d92 [ ' d9i f[ f{9,) smiq,9,) 
Jo Jo Jo 

.(11) 

vanish for n odd and \qj\ < Pj + 1. These coefficients 
are split up further by inserting the Fourier series (jlOp 
for f{9) consisting of terms sin(r9) (r odd multiple of 
+ 1). Twice the product of two sine functions is the 
difference of two cosines whose arguments are sum and 
difference of the sine arguments. Thus we want to show 



0= / d9n. 

Jo JQ Jo 



6*2 " 



d92 I d9iY[cos{{r,+qj)9,) (12) 



where rj is an odd multiple of iV -I- 1 and 

n n 

E \li I < Efe + 1) = " + ^" ^ ^- (13) 

j=o j=0 

It is easy to perform the first two integrations in (I12p 

Q 2 

cos((r3+<73)^j) [ ' d92 [ ' TT cos((rj +g,)0,) (14) 
Jo Jo 

analytically yielding a lengthy sum over terms of the form 
cos((r3 + q'^)93) where rj, is still an odd multiple of A^-|- 1 
and \q'^\ < \qi \ + \q2\ + \q3\. Hence the structure of the 
expression on the r.h.s. of (I12p is preserved, but n is 
lowered by two. This procedure is iterated till n — 1 and 
we arrive at 



0== / d9 cos{{R + Q)9) 
Jo 



(15) 
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because R is an odd multiple of + 1 and \Q\ < N, thus 
|i? + Ql G N. This concludes the derivation. 



III. LONGITUDINAL RELAXATION 

The above derivation holds also for the odd powers of 
longitudinal relaxation as observed for constant Hamil- 
tonians beforei^. The Hamiltonian studied is H{t) = 
Do(i) + Di{t) where 

Do{t) := Hk{t)+<jMt) (16a) 
Di{t) = ■ A^{t) = cr^A^it) + ayAy{t). (16b) 

The TT pulses are applied around the spin z axis so 
that the switching function appears for Di, i.e., H{t) ~ 
Do(i) + Di{t)F{t). Then all the above steps for dephas- 
ing can be repeated identically on substituting Dq 
and azAz — >■ Di. Thus we know that the UDD sequence 
of N pulses makes all odd powers up to N in Di vanish. 
Thus it efficiently suppresses longitudinal relaxation also 
for time dependent Hamiltonians. 



IV. APPLICATIONS 

The first application, of course, is the use of the UDD 
for Hamiltonians in a certain reference frames which im- 
ply that the Hamiltonians are effective ones with some 
time dependence. Examples are a rotating frame, in 
which a magnetic field is not fully compensated, or an 
effective Hamiltonian in which fast modes have been av- 
eraged by the help of Magnus expansions or they are 
treated in an interaction picture. In the cases, where 
the resulting time dependence can be considered to be 
sufficiently smooth (which need not be always true) the 
previous results establish the applicability of the opti- 
mized sequence UDD, in spite of the time dependence of 
the Hamiltonian. 

The second application concerns dynamic decoupling 
for general decoherence. The nesting of pulse sequences 
about perpendicular spin axes makes it possible to elimi- 
nate all possible couplings between a qubit and its en- 
vironment provided the expansion in time is possible, 
see Fig. [T] A first proposal used iterative concatenation 
(CDD) without optimization leading to a fast growing 
number of pulses proportional to 4^ if terms up to 
should be eliminated^. If one uses concatenation on the 
secondary level, but optimized UDD on the primary level 
(i.e. CUDD) the number of pulses grows only like 2^ im- 
proving by a square rooliii. On the primary level, CUDD 
suppresses longitudinal relaxation by tt pulses about 
the z axis such that the time evolution due to a Hamil- 
tonian H ~ Dq + Di (for Di see ([H])) without explicit 
time dependence is reduced to the time evolution due to 
an effective Hamiltonian 



ff<=«(t) = Df{t) 



(17) 



up to correction of the order 0(Tj^^+^). For a discus- 
sion of the smoothness of H°^(t) we refer the reader to 
the Appendix. The time dependence of H'^^{t) was seen 
previously as the decisive obstacle to apply an optimized 
UDD sequence again on the secondary level^^. Below we 
show explicitly that H°^{t) is indeed time dependent. 

Very recent numerical data^^, however, indicates that 
the use of a UDD sequence also on the secondary level 
with pulses about zj_ is in fact a very efficient way to sup- 
press general decoherence. This nesting of two perpendic- 
ular UDD sequences of tt pulses about z for each of the 
iVj_ + 1 intervals of a UDD sequence of N±_ vr pulses about 
a perpendicular axis is called quadratic dynamic decou- 
pling (QDD)i?. highlighting the case = Nj_. This 
choice is advantageous if longitudinal relaxation and de- 
phasing are of similar magnitude. 

The above derivation of the UDD as optimized se- 
quence for dephasing and longitudinal relaxation for time 
dependent Hamiltonians provides the analytic foundation 
of the QDD proposed by West et alJ^. In spite of the 
time dependence that the effective Hamiltonian H'^^{t) 
acquires under the primary UDD (suppressing longitu- 
dinal relaxation up to T^'+^) the secondary UDD (sup- 
pressing dephasing up to T^^+^) will still work to the de- 
sired order given by the number of pulses^. The reason 
is that an analytic time dependence, quite surprisingly, 
does not spoil the analytic properties of the optimized 
UDD sequences. 

I I l|l I I I I l|l I I >■ CUDD 

; AT time 

secondary level^ 

I I l| I L 1 |l I I > QDD 

Tlx Ttime 



FIG. 1: (Color online) Examples of nested DD suppressing 
general decoherence. Here the primary tt pulses rotate about 
the z axis and the secondary ones about x. But any pair of 
perpendicular axes will dc'^^. In practice, the primary level 
should compensate the faster decoherence. The duration of 
the primary sequences is Tp, the total duration T. 



V. TIME DEPENDENCE AFTER THE 
PRIMARY DD 

Applying tt pulses about the spin axis z to H = 
Do + Di without explicit time dependence converts 
H to H{t) = Do + DiF{t) where F{t) = ±1 is 
switching sign at the instants of the pulses^^, see 
Fig. [T] To find the time dependence of the time 
evolution t/(Tp) due to H(t) the Magnus expansion 
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f/(Tp) = exp{-irp(if«+ff(2)+^(3)) +0{T^)} is 
usediii^. The terms are powers in Tp] they read TpH^^^ — 

and 

TpH<'^^= \ r r r {[Hih),[H{h),H{t,)]] 
Jo Jo Jo 
+ [H{t3),[H{t2),H{ti)]]}dt,dt2dh. (18) 

The first order term iJ*^^^ contains the integral Ii := 
Jo^ F{t)dt which vanishes for any reasonable DD se- 
quence with at least one pulse, thus one has H^^^^ = Dq 
without time dependence. 

The second order term H^^^ involves the commutator 
[H{ti),H{t2)]. We find 

rpi/(2) ^ 1^(2) (^(i^) _ ^(i^)) diidt2, (19) 

^ Jo Jo 

where the operator t)*^^^ induces dephasing 

ry('' := J2 ^^{[Hb,A,]^^e^'^{A,,A,}), (20) 

i,j=x,y 

with e*-'^ being the Levi-Civita operator and {, } the anti- 
commutator. Since the r.h.s. of (IT^ is linear in F{t) we 
know from the general derivations given above or in Ref. 
[IB that the double integral vanishes if a UDD sequence 
with two or more pulses is applied since the total order 
of the term is T^SS. 

The third order finally yields a non-vanishing time de- 
pendence. There is a linear and a quadratic term in F{t) 



For completeness we compute the corresponding oper- 
ator [Z3i, 77*^^-'] which introduces corrections to dephasing 
and to the bath dynamics as expected 

+ io,_ {A„e'-'"^[i7b,A,]+z{A,,AJ}). (22) 

Obviously, these terms do not vanish except for very spe- 
cial choices of bath i/b and coupling operator A. This 
completes the derivation that the effective Hamiltonians 
are indeed time dependent. Thus the analytic argument 
in Ref. fisl for the validity of the QDD does not hold. 

VI. SUMMARY 

We extended the derivation of the optimized proper- 
ties of Uhrig dynamic decoupling to initial Hamiltonians 
with time dependence. First, this establishes the applica- 
bility of optimized dynamic decoupling also for the large 
class of effective Hamiltonians which inherit an explicit 
time dependence from a special reference frames or from 
the treatment of fast modes in the interaction picture, 
by average Hamiltonian theory, or by integrating them 
out. Second, our finding provides the analytic reason for 
the advantageous properties of quadratic dynamic decou- 
pling for the suppression of general decoherence including 
both dephasing and longitudinal relaxation. This scheme 
was recently proposed by West et al^. Thus the road 
is paved for a much broader applicability of optimized 
dynamic decoupling. 



^[h.j[Do,v^'^]+h.2[Di,f,^'^]} (21a) 



h,i r {F{h) + F{h) -2F{t2))dhdt2dh 

Jo Jo Jo 

(21b) 

ha - {2F{t,)F{t3) 
Jo Jo Jo 

- F{t2)F{ti)+F{t2)F{t3)))dhdt2dt3. (21c) 



We do not give the first commutator in (|21ap explicitly 
because /a^i is linear in F{t). Hence it is zero for any 
UDD with three or more pulses. But 1^,2 is quadratic 
in F{t) so that the UDD sequence does not make any 
statement on its value. Hence it will generally be finite. 
Indeed, we verify numerically for Nz pulses the examples 



I^--^ = -0.03033r3 and I^i''^ = -0.00668T3 for 



3,2 



-0.00884r3 and /g^J"^ = -0.00524T3 



odd and !■ 

for Nz even. This implies a finite quadratic dependence 
of the effective Hamiltonian H'^^ which we aimed to es- 
tablish. 
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Appendix A: Effective Hamiltonian from the 
primary level of QDD 

Here we discuss in more detail how H°^{t) arises from 
the time evolution on the primary level in QDD. In par- 
ticular, we discuss why this operator is smooth or even 
analytic in t although it arises from the integration of 
H{t) = Dq + DiF{t) where F{t) is a switching function 
which changes sign at the pulses of the primary level. 

Let us denote by U{t) the time evolution operator on 
the primary level, i.e., within one of the intervals of the 
secondary level, see Fig. [T] It reads 



U{t) = Texp {-if H{t')dt') 
Jo 



(Al) 



and the domain of integration is illustrated in the upper 
panel of Fig.[2]by the range of the switching function F{t) 
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FIG. 2: (Color online) Ranges of integration (range of solid 
curve) and behavior of the switching function F{t) for the 
calculation of U{t) (upper panel) and of the effective f/°^(t) 
(lower panel). 

encountered. If we define the dimensionless relative time 
t :— t/Tp and the corresponding Hamiltonian 
H(Tpt) the previous equation becomes 

/•*/^p 

U{t) = Texp { -iTp H'''\i)dt) (A2) 

where T stands for the time ordering according to the 
relative time t. Eq. (jA2p is given here for comparison 



with the subsequent effective time evolution operator. 

To consider what happens on the secondary level, we 
need W'^it). Its defining property is W'^iTp) := UiTp). 
But if its argument is varied Tp ^ t it is implied that 
the primary sequence is scaled accordingly. That is its 
decisive difference to U{t) which is illustrated in the lower 
panel of Fig. [2l Note that this scaling is exactly what is 
done when the primary sequences are applied in each of 
the time intervals of varying duration of the secondary 
level, see Fig. [1] 

In a formula, J7°* {t) is given by 

U''^{t)=fexp{-tt[ H"'Hi)dt). (A3) 
Jo 

From this equation it is obvious that (t) is generically 
smooth in the variable t. If the parts of the Hamiltonian 
are bounded U'^^{t) is analytic because the exponential 
is analytic. 

The final step is to define the corresponding effective 
Hamiltonian H'^^{t). As usual the Hamiltonian is re- 
trieved from the time evolution as its infinitesimal gen- 
erator 

W^{t) i [dtW^it)] [U''^[t)\ ^ . (A4) 

If C/°*^(i) is analytic, then H^'^it) is it as weU. If C/°*^(t) 
can be expanded up to and including the power t^"^^, 
then i/°*^(t) can be expanded up to and including the 

4.N 

power t . 
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